function eqinox

global imode lmode

lmode = imode;

if (imode <= 1)

    imode = imode + 2;

end